Abstract 

The problem of radio wave reflection from an optically thick plane 
monotonous layer of magnetized plasma is considered at present work. 
The plasma electron density irregularities are described by spatial 
spectrum of an arbitrary form. The small-angle scattering approx- 
imation in the invariant ray coordinates is suggested for analytical 
investigation of the radiation transfer equation. The approximated 
solution describing spatial-and-angular distribution of radiation re- 
flected from a plasma layer has been obtained. The obtained solution 
has been investigated numerically for the case of the ionospheric ra- 
dio wave propagation. Two effects are the consequence of multiple 
scattering: change of the reflected signal intensity and anomalous re- 
fraction. 
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1 Introduction 



Basic goal of the present work consists in derivation and analysis of the 
transfer equation solution describing spatial-and-angular distribution of radio 
radiation reflected from a plane stratified layer of magnetized plasma with 
random irregularities. 

The radiation transfer equation (RTE) in a randomly irregular magne- 
tized plasma was obtained in the work |l[ under rather general initial assump- 
tions. In particular, the medium average properties were assumed smoothly 
varying both in space and in time. In the work the radiation energy 
balance (REB) equation describing radiation transfer in a plane stratified 
layer of stationary plasma with random irregularities has been deduced. The 
invariant ray coordinates, allowing one to take into account, by a natural 
way, refraction of waves and to represent the equation in the most simple 
form, were used there. In the work || it was shown that the REB equa- 
tion is a particular case of the radiation transfer equation obtained in [|TJ 
and can be deduced from the latter by means of transition to the invariant 
ray coordinates. REB equation, thus, allows one to investigate influence of 
multiple scattering in a plane stratified plasma layer on the characteristics of 
radiation. In particular, it enables one to determine the spatial-and-angular 
distribution of radiation leaving the layer if the source directivity diagram 
and irregularity spatial spectrum are known. A few effects which require 
of wave amplitudes coherent summation for their description (for example, 
phenomenon of enhanced backscattering) are excluded from consideration. 
However, the multiple scattering effects are much strong cr, £lS 8b rule. This is 
true, in particular, for the ionospheric radio propagation. 

The numerical methods of the transfer equation solving developed in the 
theory of neutron transfer and in the atmospheric optics appear useless for 
the REB equation analysis. They are adapted, basically, to the solution of 
one-dimensional problems with isotropic scattering and plane incident wave. 
In of magnetized plasma the presence of regular refraction, aspect- 

sensitive character of scattering on anisometric irregularities and high di- 
mension of the problem (the REB equation contains two angular and two 
spatial coordinates as independent variables) complicate construction of the 
effective numerical algorithm for its solving. In this situation it is expedient 
to solve the REB equation in two stages. The first stage consists of obtain- 
ing of the approximated analytical solution allowing one to carry out the 
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qualitative analysis of its properties and to reveal of its peculiarities. At 
the second stage the numerical estimation methods can be applied to the 
obtained analytical solution. This approach has been realized at our work. 

We begin the present paper from a detailed exposition of the invariant 
ray coordinates concept. Then possibility to use of the small- angle scattering 
in the invariant coordinates approximation is discussed. Two modifications 
of the REB equation solution are obtained. The analysis of the obtained 
solutions (both analitical and numerical) concludes the paper. 

2 Invariant ray coordinates and the radiation 
energy balance equation 

It is convenient to display graphically the electromagnetic wave propagation 
in a plane-stratified plasma layer with the aid of the Poeverlein construction 
H, We shall briefly describe it. Let the Cartesian system of coordinates 
has axis z perpendicular and the plane xOy parallel to the plasma layer. We 
shall name such coordinate system "vertical" . It is assumed that the vector 
of the external magnetic field H is situated in the plane zOy . Module 
of radius- vector of any point inside of the unit sphere with centrum in the 
coordinate origin corresponds to the value of refractive index rii(v , a), where 
i = 1 relates to the extraordinary wave, i = relates to the ordinary one, 
v = oj 2 /uj 2 , u 2 is the plasma frequency, u 2 is the frequency of a wave, a is 
the angle between radius-vector and magnetic field H . The refractive index 
surface corresponding to a fixed value of v and to all possible directions of 
the radius-vector represents a rotation body about an axis parallel to vector 
H (see fig. 1). 

Convenience of the described construction (in fact, this is an example 
of coordinate system in space of wave vectors k) is become evident when 
drawing the wave trajectory: it is represented by a straight line, parallel to 
the axis z. This is a consequence of the generalized Snell law, which also 
requires of equality of the fall angle and exit angle onto/from a layer (6), and 
constantness of the wave vector azimuth angle (ip). Note, that the crossing 
point of a wave trajectory with the refractive index surface under given value 
of v determines current direction of the wave vector in a layer (it is anti- 
parallel to a radius-vector) and current direction of the group speed vector 
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(it coincides with the normal to the refractive index surface). The projection 
of a wave trajectory onto the plane xOy is a point which radius-vector has 
module sin# and its angle with relation to axis x equals to tp. Thus, the 
coordinates 9 and tp define completely the whole ray trajectory shape in a 
plane layer and outside of it and are, in this sense, invariant on this trajectory. 

Radiation of an arbitrary point source of electromagnetic waves within 
the solid angle 9 -j- 9 + d9; tp -v- <f + dtf corresponds to the energy flow in the 
fc-space inside of a cylindrical ray tube parallel to axis z with cross section 
sm9d(sin9)dip = sin9 cos9d9dp. In case of regular (without random irregu- 
larities) plasma layer this energy flow is conserved and completely determined 
by the source directivity diagram: 

P(z,p,9,p)=P (p,9,p) (1) 

where P is energy flow density in the direction determined by angles 9, p 
through the point p on some base plane situated outside of the layer parallel 
to it (in the ionosphere case it is convenient to choose the Earth' surface as 
the base plane), z is distance from the base plane (height in the ionosphere 
case). We shall assume in the present paper that function z(v ) is monotonous 
in the region of wave propagation and reflection. If random irregularities are 
absent and source of radiation is point, variable p in (|IJ) is superfluous, 
as the matter of fact, since unequivocal relation between it and angles of 
arrival of a ray 9, p exists. When scattering is present the radiation energy 
redistributes over angular variables 9, p and in space what is described by 
variable p. The value of P satisfies in this case to the equation of radiation 
energy balance B |3| : 



f z P(z, 9, p, p) = J{-P(z, p, 9, tp) sin 9 cos 9C~\z- y 9, tp). 
■a [« (9, <p) , fa (9, <f) ; a (0', if') , (3 (9', <p>)} sin a(9', p>>) 
+P 



V tV >) 



+ 



(2) 



z, p - ®(z; 9', p'- 9, (p), 9', p'\ sin 9' cos9'C- 1 {z; 9', p')- 
■a [a (9', iff) , fa (9\ iff) ;a(9,ip),(3 (9, tp)] sma(9, <p) \^\}d9'dp' 

where C(z; 9, tp) is cosine of a ray trajectory inclination angle corresponding 
to the invariant angles 9 and tp; \d(a, /3)/d(9, tp\ is Jacobean of transition 
from angular coordinates 9, p to the wave vector polar and azimuth angles 
a and j3 in the "magnetic" coordinate system (which axis Oz is parallel to 
the magnetic field); 

a [a (9, if) , Pa (9, p) ; a (9', tp') , p (9'.tp')} = a [9, tp- 9', tp') 
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is scattering differential cross section describing intensity of the scattered 
wave with wave vector coordinates a, (3 in magnetic coordinate system (cor- 
responding invariant coordinates are 9' and ip') which arises at interaction of 
the wave with wave vector coordinates «o, (3$ (invariant coordinates 9 and tp 
) with irregularities. Vector function <3>(z; 9', ip'; 9, ip) represents the displace- 
ment of the point of arrival onto the base plane of a ray which has angular 
coordinates 9' and ip' after scattering at level z with relation to the point of 
arrival of an incident ray with angular coordinates 9, ip . It is essential that 
in a plane-stratified medium the function $ is determined only by smoothed 
layer structure v(z) and does not depend on the scattering point horizontal 
coordinate and also on coordinate p of the incident and scattered rays. Note 
also that ratio 9, ip; 9', ip') = — 9', ip'; 9, ip) takes place. 

It is possible to check up that equation (2) satisfies to the energy con- 
servation law: when integrating over all possible for level z values of 9, ip 
and all p its right side turns into zero. It is natural since in absence of true 
absorption the energy inside the plasma layer does not collected. 

Analyzing expression for the scattering differential cross section in a 
magnetized plasma (see, for example, ||), it is easy to be convinced that the 
following symmetry ratio takes place: 

a [9, ip; 9', ip'} n 2 cos &' g = a [9', ip'; 9, ip] n' 2 cos $ g (3) 

where $ g is angle between the wave vector and group speed vector, n is 
refractive index. Using @ the equation (Q) can be presented as follows: 

£P(z,p,9,ip) = JQ(z;9,ip;9',ip>) 

{P(z,p- <f(z;9',ip';9,ip),9',ip')-P(z,p,9,ip)}d9'dip' 1 > 

where Q(z; 9, ip; 9', ip') = a{9, ip; 9', <//)C -1 (,2, 9, ip) sm9' \dQ! k /dQ!\ , and quan- 
tity Q (z; 9, ip; 9', ip') = Q(z; 9, ip; 9', ip') sin 9 cos 9 is symmetric with relation 
to rearrangement of shaded and not shaded variables. The equation REB in 
the form (^) has the most compact and perfect appearance. It is clear from 
physical reasons that (Q) has to have a unique solution for given initial distri- 
bution Po(p, 9, ip). The obtained equation can be directly used for numerical 
calculation of the signal strength spatial distribution in presence of scatter- 
ing. However, as it was noted at introduction already, this approach leads to 
essential difficulties. Subsequent sections describe the method of construction 
of the energy balance equation approximated analytical solution. 
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3 Small-angle scattering approximation in the 
invariant ray coordinates 

Let us consider the auxiliary equation of the following kind, which differs 
from (4) only by absence of the dash over variable uj marked by arrow: 



dz 



-P{z,p,») = jQ(z;^'){p(z,p + ^J), 1 )-P(z, M }^ 

(5) 

where designation uj = {9, ip} , du = dOdip has been used for the sake of 
compactness. Equation (||) can be easily solved analytically by means of 
Fourier transformation over variable p. The solution has the following form: 

P (z,q,u) = P (q,uj)S(z,0;q,u) (6) 

where Po(q,uj) is the Fourier image of the energy flow density of radiation 
passing the layer in absence of scattering and the value of S is defined by the 
expression 

S(z 2 , Zl , q, uj) = exp { / dz' f du'Q(z'; u; uj') \ e ^ z '^ - ll 1 (7) 



One should note that integration over z in this and subsequent formulae, in 
fact, corresponds to integration along the ray trajectory with parameters 9, if. 
The area of integration over u' includes rays which reflection level h r {uj') > z. 
Let us transform now equation (^) by the following way: 

£P(z, p,u)=J du'Q(z; u; uj') {p(z, p+ <S(z; uj; uj'),u) - P(z, p,uj)} + 
+ J duj'Q(z;uj;uj') {P(z, p+ §(z;uj;uj'),uj') - P(z, p+ $(z; uj; uj'), uj)\ 

(8) 

Its solution will be searched for in the form 

P{z,p,u)=P {z,p,uj)+X(z,p,uj) (9) 

Thus, auxiliary equation (H) allows to present the solution of the equation 
(H) in the form (^). This is an exact representation while some approximated 

expressions for quantities P and X are not used. 
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By substituting of (|9|) into the equation (f|) one can obtain the following 
equation for the unknown function X : 



dz 



X(z,p,uj) — f du'Q(z;uj;uj , ){p z,p + $(z;uj;uj'),u}' 



P 



z,p + &(z;lu;uj'),uj }+ 



(10) 



+ J du'Q(z;u;u'){X z,p + $(z; o>; u/), u/ — X(z, p, a>)} 



We shall assume now that the most probable distinction of angles u' and to 
is small. The heuristic basis for this assumption is given by analysis of the 
Poeverlein construction (fig. 1). It is easy to be convinced examining the 
Poeverlein construction that scattering near the reflection level even for large 
angles in the wave vector space entails small changes of the invariant angles 
8, (p. This is especially true for irregularities strongly stretched along the 
magnetic field (in this case the edges of scattered waves wave vectors form 
circles shown in fig. 1 as patterns A and B). One should note also that the 
changes of invariant angles 9, if are certainly small if scattering with small 
change of a wave vector direction takes place. This situation is typical for 
irregularity spectra, in which irregularities with scales more than sounding 
wave length dominate. Thus, the small-angle scattering approximation in the 
invariant coordinates has wider applicability area than common small-angle 
scattering approximation. 



Scattering with small changes of 8, ip entails small value of 
bllows directly both from sense of this quantity and from the 



That 



'act what 



$(z,uj,U}) = 0. Let us make use of that to carry out expansion of quantity 
X at the right side of the equation fllQ|) into the Taylor series with small 
quantities &' — uj and $ . Note that making similar expansion of function 
P at the initial equation (Q) would be incorrect since function P may not 
to have property of continuity. For example, in case of a point source, Pq 
is a combination of (^-functions. As it will be shown later, the function 
X is expressed through P by means of repeated integration and, hence, 
differentiability condition fulfils much easier for it. 

Leaving after expansion only small quantities of the first order, we obtain 
the following equation in partial derivatives: 



d_ 

dz 



d 

X(z,p,uj) - A UJ (z,u)—X(z,p,u) + A t 



d_ 

dp 



X{z,p,uj) 



/( z,p,u) 
(11) 
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where 



A uj (z,uj) = j duj'Q{z;uj;uj l ){uj' — u)\ 

Ap{z, uj) = J cIuj'Q(z;uj;uj')$(z,uj,uj'); 
/( z, P» = I duj'Q(z; uj] uj')- 



P 



z,p + §>(z;uj;uj'),uj' — P z,p + $(z;uj;uj'),uj 



Here is the characteristic system for the equation (p"T|): 

dX _ dp . , .did 

~T = /( z,p,u)] — = Ap{z,u))] — = -A UJ (z,u) 
dz dz dz 

and initial conditions for it at z — 0: 

X = 0; () = p; uj = ujo. 

It is necessary to emphasize the distinction between quantity p 1 , which is a 
function of z, and invariant variable p. 

Solving the characteristic system we obtain 

20 

uj = uj(z,uj ), p = p- J dz'Ap[z',uj(z',uj )] 

z 

where Zq is z coordinate of the base plane. It follows that 



-0 



X(z ,p,uj) =J dz'f{ z',p- I dz"Ap[z",uj(z",uj )},uj(z',uJo) 



(12) 



Generally, expression (|12"D gives the exact solution of the equation fliTf). How- 
ever, since we are already within the framework of the invariant coordinate 
small-angle scattering approximation which assumes small value of A w (z, uj), 
it is possible to simplify the problem a little. Assuming A^ = and omitting 
index at invariant coordinates uj, we are earning to the following approxi- 
mate representation for function X : 



20 



X(z ,p,uj) = / dz'{P z',p+ $(z';uj;uj') + D(z ,z',uj),uj' 



P 



(13) 



z',p+ $(z';uj;uj') + D(z , z',uj),uj } 
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where D(z%, Zi, lu) 



dz' J duu'Q{z'\ lu; Lu')$(z r , lu, lu'). 



Thus, in the invariant coordinate small-angle scattering approximation 
the solution of the REB equation (|J) is represented as a sum of two terms 
(see (|^)), the first of which is 



P (zo,P» 



(27T)' 



Jd 2 qP (q,Lu) 



exp < iqp+ / dz' f du'Q(z'; lu; lu') 



(14) 



where 



/ d 2 q Po(q,w) exp(iqp) = P (p,Lu), and the second one is given 



by expression (|T3|). 

The solution can be presented in the most simple form if one uses again 
the smallness of quantity $ and expands the second exponent in the formula 
( |T4] ) into a series. Leaving after expansion only small quantities of the first 
order, one can obtain: 



2(1 



P(z ,p,u) £ P p + D{z ,0,u),uj + / dz' J du'Q(z';u;u') 



■{P z',p+ <f>(z';u;uj') + D(z ,z',u) + D(z',0,u'),u' 



(15) 



P 



z',p+ $(z';uj;uj') + D(z o ,0,uj),uj } . 



The last operation is the more precise the faster value of Po(q, to) decreases 
under |g| — > oo. The solution of the radiation energy balance equation 
obtained in the present section in the form @, (0), (|i~3l), or in the form (|T5|), 
expresses the spatial-and-angular distribution of radiation intensity passing 
layer of plasma with scattering through the spatial-and-angular distribution 
of the incident radiation, that is, in essence, through the source directivity 
diagram. 



4 Alternative approach in solving the REB 
equation 

The REB equation solving method stated in the previous section is based on 
representation of quantity P(z, p, to) as a sum of the singular part P (z, p, uj) 
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and the regular one X{zq, p,u>). Regularity of the X(zq, p, uj) has allowed 
one to use the expansion into the Taylor series over variables p and uj at the 
equation ( [TO| ) right side and to transform the integral-differential equation 
fllOD into the first order partial derivative differential equation (pTT|) . 

However, the stated approach is not the only possible. The REB equation 
can be transformed right away using Fourier-representation of the function 
P(z,p,u): 



P(z,p,uj) 



(2tt)' 



d 2 qP{z, q, uj) exp{iqp) 



(16) 



Substitution of (|T^) into (|) gives the following equation for quantity P(z, p, uj): 

—P{z,q,oj) = J dw'Q{z; oj; oj') {P{z, q, oj') exp iq${z; uj; uj') — P{z,q,oj)} 

(17) 

The quantity P(z, q, oj) is a differentiable function even when P{z, p, uj) has 
peculiarities. Therefore, in the invariant coordinate small-angle scattering 
approximation it is possible to use the following expansion: 



P{z, q, uj') - P{z, q, uj) + dP{Z d J ,U}) (u/ - uj) 



(18) 



Substituting (|TB|) in (O) we obtain the partial derivative differential equation 

d ~ d 

—P{z,q,oj)- A {z,q,oj)—P{z,q,oj) - P{z,q,u) S {z,q,u) = (19) 
where 



S {z,q,oj) = / duj'Q{z;uj;uj') 



,JqQ(z;u>;uj') j 

A{z,q,u>)= J doj'Q{z;oj;oj') exp iq&{z; oj; oj') {uj' — uj) . 

The characteristic system 

duj ~ „ ^ . dP 



dz 



dz 



~-S {z,q,uj)P{z 1 q,uj) 



(20) 



with initial conditions P = Po{q,uj), uj = ujq at z = has the following 
solution: 



P{z, p, uj) = ——^ J d 2 qP {q, uj ) exp j iqp+ J dz' S [z f , q, uj{z', q, 



(21) 
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This solution of the REB equation turns into the expression ( |T4] ) for P when 

A (z, q, u) — ► . But the latter limit transition corresponds to the invariant 
coordinate small-angle scattering approximation used in the previous section 
under derivation of ([13|) and subsequent expressions. Let us note, however, 
that in (pT|) , in contrast with @, any additional terms do not appear. It 
allows one to assume that in used approximation the ratio 

X(z,p,u) < P(z,p,u) (22) 

is fulfiled. Additional arguments to the benefit of this assumption will be 
presented in the following section. 



5 Analysis of the solution of the REB equa- 
tion 

We shall show, first of all, that the obtained solution satisfies to the energy 
conservation law. For this purpose it is necessary to carry out integration 



of the left and right sides of (15) over lu and p multiplied them previously 
by sin 9 cos 9. The area of integration over angles is defined by the condition 
that both wave uj and wave lu' achieve the same level z (since at level 
z their mutual scattering occurs). To satisfy this condition one should add 
factors O [h r (u) — z] and [h r (uj') — z] to the integrand expression, where 
Q(x) is the Heviside step function, h r (uj) is the maximum height which can 
be reached by a ray with parameters 9, if. Now integration can be expanded 
over all possible values of angles, i.e., over interval -j- 7r/2 for 9 and over 
interval 4- 2tt for ip. Then, (JIB]) becomes 



/ P(ui) sin 9 cos 9du = J Pq(oj) sin 9 cos 9duj+ 

z 

+ / dz' JduJ duj'Q [h r (u) - z'\ [h r {u') - z'\ Q (z'; a, <J) [P (u') - P (u)} 



(23) 

where P(uj), P (u) is a result of integration of P(z ,p,u) and P (p,uj) 
correspondingly over variable p. 

Due to antisymmetry of the integrand expression with relation to rear- 
rangement of shaded and not shaded variables, the last term in (|23|) is equal 
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to zero. Thus, equation (|23| ) reduces to 

J P(z , p, to) sin9 cos9dtud 2 p = J Po(p, oj) sin9 cos9du>d 2 p (24) 

expressing the energy conservation law: the radiation energy full flow through 
the base plane remains constant regardless of scattering, as it should be in 
case of real (dissipative) absorption absence. It is not difficult to check that 
parity (|24] ) is valid for the exact solution in the form (||) and also for the 
solution in the form (pl|). 

With relation to the solution in the form (H) the carried out discussion 
discovers one curious peculiarity. It appears that the radiation energy com- 
plete flow through the base plane is determined by the first term ^P^j . The 

second one (X) gives zero contribution to the complete energy flow. 

Let us investigate in more detail the structure of quantity X(z,p,u>) in 
the invariant coordinate small-angle scattering approximation. Proceeding 
to the Fourier-representation in the expression (|13|) produces 



X(zq, q,u>) —J dz' J duj'Q(z'; u, u>') 
^ o 

•expjzg <£>(V; lu; lu') + D(z , 0, lu) j 



P {z',q,u')- P (z',q,u) 



Employing regularity of function P (z,q,u), the last expression can be 
written as 



X(z ,q,u) =J dz 
o 



,d P (z',q,u) 



~ r — * 

A (z , q, lu) exp iqD(z , z , to 



where quantity A (z,q,cu) is defined by (|19|) . Thus, it becomes evident 

that limit transition A (z,q,u) — > entails also X(zq,p, uj) — > . This 
property has been established in previous section with the aid of comparison 
of two variants of the REB equation solution. Now we can see that its 
presence is determined by structure of quantity X(z , p, to). 

Results of the present section give the weighty ground to believe that 
the radiation spatial-and-angular distribution is determined basically by the 
first term in the solution (|9|). The second term represents the amendment to 
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the solution which can be neglected in the invariant coordinate small-angle 
scattering approximation. This statement validity can be checked under 
detailed research of properties of the obtained REB equation approximated 
solutions by numerical methods. 



6 



Technique of numerical calculation of mul- 
tiple scattering effects 



We shall proceed from the obtained solution ( |15|) where only the first term 
has been retained. In the considered approximation the multiple scatter- 
ing results in deformation of the radiation field reflected by a plasma layer 
without change of a kind of function describing intensity spatial and angular 
distribution. If only single ray with parameters 9q (p), <po ifi) comes into each 
point p onto the base plane when reflecting from a regular plasma layer (it 
will be so, if the source is dot-like and the frequency is less than critical one 
for this layer), for function P it is possible to use expression of a kind 



P (p, 9, if) =P (p) 5 [cos 9 - cos 9 (p)] 5{ip-p (p)] 



(25) 



where quantity Pq (p) has the meaning of energy flow at the point p in 
absence of scattering. Substituting of (|25| ) into ( |i~5D and making integration 
over angles one can obtain for the energy flow at the point p 



P (P) =Po 



9, 



Op 



-9 



Op 



P + D{9 uVl ) 



Pop 
D, 



p + D{6i,<pi 
p + D{9 l) ^ l 

, (01, (fx) (f p 



D 6 (9 1 ,< fl )}- 
D v (9 1 , Vl )}- 
p + D{9 1 ^ 1 )} D 6 



(26) 



h,<Pi) 



where D (9, if) = D (z , 0; u>), subscripts p, 9, ip mean derivatives with corre- 
sponding variables, and 9\, (fx represent new arrival angles of a ray. 

Expression ( p6|) uses explicit dependencies of arrival angles #o, y?o of a 
ray reflected from a regular plasma layer on position p which are usually 
unknown. As a rule, the dependence of coordinates p on 9q and <po can be 
expressed in an explicit form: p = po (6*o, v^o) , where po (9, ip) - point of arrival 
onto the base plane of a ray with invariant angles 9 and ip. Expressing #op 
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and ifQp in (|26| ) via — and — , we obtain new representation for quantity 
P (P)-- 



d(po x ,POy) 



d (pox — D x , po y — Dy 



(27) 

New arrival angles #1 and tpi of a ray can be found by solving the algebraic 
equation system 

p-p (e 1 ,<p 1 ) + D(e 1 ,<p 1 ) = o (28) 



According to (26) , an observer being at point p discovers two effects 
connected with scattering in a plasma layer: change of the wave arrival angles 
and change of the received signal intensity. The analytical results of the 
present work are valid in a common case of magnetized plasma. However, we 
shall consider below the case of an isotropic plasma linear layer, as it displays 
the task basic features at relative simplicity of numerical calculations. The 
simplification is mainly due to possibility to define the ray trajectories in 
the analytical form. For a numerical estimation of the effects described by 
the main term of the REB equation solution it is necessary, first of all, to 
concretize the kind of the AN/N irregularity spectrum. We shall choose the 
spectrum of the following kind: 



F{R) = C a {i + k\/4 a ) v '\(k 



(29) 



where k± and Ku are vector k (irregularity spatial harmonic) components 
orthogonal and parallel correspondingly to the magnetic field force lines, 
K 0J _ = 2ir/l ±, l ± is the spectrum external scale, 8 (x) is the Dirac delta- 
function, 



C A = 5 



R' 



r fr/2) 

27TKq_ | _ 



r 



V-2 N 



("~2)/2 



K( u -2)/2 (Rk>0±) 



is a normalization constant, T (x) is the Gamma-function, Kp (z) is the Mac- 
donald function [^]. The 5r quantity characterizes the level of irregulari- 
ties AN/N. In mathematical theory of random fields it corresponds to the 
structural function of the irregularity field for the scale length R ||. The 
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considered model of a power-like small-scale irregularity spectrum is used 
in many areas of modern physics. Both at the ionospheric F region and at 
the tokamak plasma the irregularities are strongly stretched along the mag- 
netic field. We preserve this feature of irregularities even when dealing with 
isotropic plasma model. Let the plasma layer with a linear dependence of 
electron density on depth (dz/dv = H) be located at distance h from the 
radiation source. The scattering cross-section for isotropic plasma is ||: 

a = ^kyF(AK) (30) 

where ko = u/c, F (k) is the irregularity spatial spectrum and An is the scat- 
tering vector. For the case of infinitely stretched irregularities the scattering 
vector longitudinal and transversal components is defined by the expressions 

A/cm = fcon (cos a — cos a') ; 

Ak\ = k 2 n 2 (sin a cos (3 — sin a' cos (3') 2 + (sin a sin (3 — sin a' sin f3') 2 

where notation of the first section is used and isotropic plasma refractive 
index is n = a/1 — v. In a linear layer of isotropic plasma the vector $ 
components can be presented as 

$ x 0; Ok, Vk, 9' k , <p' k ) = f (e' k ) cosp'k - f (0 fc ) cos^ fc ; 

<t> y (v; 6 k , (p k , 9' k , (p k ) = f (6' k ) sin <f' k -f {B k ) sin (p k ; 

where 



/ (9k) = 2Hn sm8 k (yl — n 2 sin 9 k + n cos 0^1 + ho sin 9 k /y 1 — n 2 sin 9 k , 

angles 9 k , ip k , 9 k , <p' k are current polar and azimuth angles of the wave vectors 
of incident and scattered waves correspondingly at the "vertical" coordinate 
system. 

Because of 5-function presence in the irregularity spectrum (^9|), the nu- 
merical estimation of expression for D (9, (p) is reduced to calculation of dou- 
ble integral over v and one of angles. It is convenient to proceed to integration 
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over angle P' at the "magnetic" coordinate system. As result we obtain 



cos 2 9 



D x = yk 3 H 2 C A 



dvv 2 



1 



2- 



cos 2 9 



d(3{{smai cos/?- 



he 



2 (y/v + Wl (l + v) + WtVT^) + ~§( v + W i (! - v )) 

ho sin 9 

H y/1 — V COS 9 



-1/2 



-COS (f 



sin 



(cos# ± Vcos 2 9 — v ) + 



2.-2 ^ 



f ) sin 2 «i sin 2 



-i//2 



(31) 



+ 



+ similar term in which aci, Pihave been replaced by a 2 , P2} 

where W\ = cos aci cos 7 — sin a x sin /3 sin 7, angles and a 2 ,/? 2 are po- 

lar and azimuth angles of the incident wave wave vector at the "magnetic" 
coordinate system for the trajectory ascending and descending branches cor- 
respondingly. These are connected to the invariant angles by the relations 



cosai >2 = ± cos 71 



'cos 2 6 — v sin 9 



1-v 



sin (p sin 7; 



sin Pi 2 sincti 2 



sma 12 = 
sin 9 



cos /?i j2 sina li2 



1 — cos 2 «i j2 ; 
sin ip cos 7 ± 



cos 2 9 — v 
1-v 



■ sm7; 



sin 9 



COS(p 



s/T-V 

where the top signs correspond to the subscript 1, the bottom ones corre- 
spond to the subscript 2. Expression for D y is derived from the expression for 
D x by replacement of underlined factors sin qji )2 cos P with (cos a ii2 sin 7 + sin a ii2 sin P cos 7) 
and cos if with simp correspondingly. Two terms in fl3T|) correspond to the 
trajectory ascending and descending branches of a ray with coordinates 9, ip. 
The equation system fl28|) numerical solving was performed using the New- 
ton's globally converging method described, for example, in M. 
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7 Calculation results for isotropic plasma 



The calculations were carried out for the following set of parameters: ho = 150 
km, H = 100 km, v = 2.5, Iq± = 10 km, R = 1 km, Sr = 3- 10~ 3 (ionospheric 
irregularity level characteristic of night quiet conditions), frequency / = 5 
MHz, the angle of irregularity inclination ("magnetic field" inclination) 7 = 
25°. Thus we mean the conditions of the ionosphere sounding from the 
Earth surface. The intensity attenuation and change of the arrival angles of 
a signal reflected from a layer (with relation to those values when reflecting 
from the same plasma layer without irregularities) were calculated for area 
which sizes were 800 km along the y axis (the magnetic meridian direction) 
and 400 km along the x axis. The radiation source was in the coordinate 
origin. The calculation results are represented in figures 2-4. On the 
contour map of fig. 2 the constant level lines of the received signal intensity 
attenuation (in dB) are shown. The signal attenuation calculated value is 
symmetric with relation to both the y axis (i.e. magnetic meridian plane) and 
the x axis with accuracy determined by numerical calculation errors. The 
result is not trivial: there is asymmetry in the problem conditions due to the 
irregularity inclination. The central symmetry is required by the reciprocity 
theorem [f|, |5|. Thus the obtained numerical solution is in accordance with 
the electromagnetic field general properties. This is an additional argument 
to the benefit of the conclusion about the primary significance of the first 
term in the REB equation approximate solution (|T5|). 

The main detail in fig. 2 is the region having the shape of ellipse with 
half-axes of 300 km and 60 km where significant intensity attenuation (up to 
15 dB) takes place. Outside of this region some increase of the signal intensity 
(in comparison with its value in absence of irregularities) is observed. This 
is quite natural result because at an nonabsorbing medium the complete 
radiation energy flow is conserved and scattering results only in its spatial 
redistribution. At larger distance from the source the intensity change aspires 
to zero. The transition from the region of the reflected signal attenuation 
to the region of amplification is of sharp character. That is, probably, a 
consequence of approximation used under transformation of the expression 
(|I4]) to the expression (|15j) fist term. One can expect that retaining of the 
higher order terms in the series expansion will result in smoothing of the 
above transition. 

The second of multiple scattering effects (i.e. change of the arrival angles) 
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is illustrated if figs. 3 and 4. In fig. 3 the contour lines show the absolute 
value of the polar arrival angle 9 alteration. In fig. 3 the alteration of azimuth 
angle is presented. One can see that distortion of the polar angle reaches of 
5 and maximum alteration of the azimuth angle is 90°. 

Both effects are observed in experiment. The intensity reduction of the 
vertical sounding signal reflected from the ionosphere should be interpreted 
by the observer as an additional collisionless mechanism of radio wave atten- 
uation. This phenomenon is observed under natural conditions bringing the 
increased values of the effective collision frequency Jl0|, [llj]. There are the 
weighty grounds to believe that it is connected to development of small-scale 
irregularities in the ionosphere. In particular, it is displayed stronger at night 
time and can reach of 10-15 dB [|TTJ . The latter figure is in accordance with 
our calculation results. The experimental data on the arrival angle change 
of a wave on a short line (when transmitter to receiver distance is about 100 
km) are reported in fT2|| . 



8 Conclusion 

In the present work the heuristic basis for use of the invariant coordinate 
small-angle scattering approximation under solving of the RTE for a magne- 
tized plasma layer is considered. Within the framework of this approximation 
two versions of the analytical solution have been obtained. They describe 
spatial-and-angular distribution of radiation reflected from a monotonous 
plasma layer with small-scale irregularities. 

The physical conclusions about influence of the multiple scattering effects 
in a layer of plasma on the spatial-and-angular characteristics of radiation are 
possible on the basis of detailed numerical research of the obtained solutions. 
Such research has been carried out in the present work for the case of isotropic 
plasma. 

It was shown that the main term of the REB equation solution for the 
radiation reflected from a plasma layer with random irregularities describes 
two effects: the signal intensity change (attenuation for the normal sounding) 
and the arrival angle change. Both effects are observed in the experiments on 
the ionosphere radio sounding. The first one is known as anomalous attenu- 
ation of radio waves at natural conditions. Note, that under the ionosphere 
heating experiments another kind of anomalous attenuation is observed: its 
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mechanism is based on the mode transformation under scattering, not on the 
multiple scattering [ 13] . The effect of arrival angle change can be interpreted 
as a mechanism of additional refraction in the ionosphere and also has exper- 
imental confirmation. These two effects numerical estimations obtained in 
the present work for parameters of a plasma layer and irregularity spectrum 
typical for mid-latitude ionosphere are in accordance with experimental data. 

The considered effects can be observed not only at the ionosphere radio 
sounding but also at sounding by electromagnetic radiation of other kinds of 
plasma with random irregularities both in natural and in laboratory condi- 
tions. 
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Figure captions 



Fig. 1. A schematic plot of Poeverlein's construction for ordinary waves. 
The refractive index surfaces for several values of v (v — uj^/uj 2 ) are shown. 
The ray trajectories in this "/c-space" are represented by straight dashed lines 
parallel to z axis. 

Fig. 2. Contour map showing the relative alteration of the intensity (in 
dB) of the signal reflected from the ionospheric plasma layer due to multiple 
scattering. Radiation source is in the coordinate origin. 

Fig. 3. Contour map showing the relative alteration of the polar arrival 
angle (in degrees) of the signal reflected from the ionospheric plasma layer 
due to multiple scattering. 

Fig. 4. Same as in fig. 3, but for the azimuth arrival angle. 
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